Measurement of local dissipation scales in turbulent pipe flow 
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Local dissipation scales are a manifestation of the intermittent small-scale nature of turbulence. 
We report the first experimental evaluation of the distribution of local dissipation scales in turbu- 
lent pipe flows for a range of Reynolds numbers, 2.4 x 10^ < Rcd < 7.0 x lO**. Our measurements 
at the nearly isotropic pipe centerline and within the anisotropic logarithmic layer show excellent 
agreement with distributions that were previously calculated from numerical simulations of homo- 
geneous isotropic box turbulence and with those predicted by theory. The reported results suggest 
a universality of the smallest-scale fluctuations around the classical Kolmogorov dissipation length. 

PACS numbers: 47.27.eb, 47.27.nf 



The dissipation of turbulent kinetic energy is driven 
by fluctuations of velocity gradients whose amplitudes 
can locally and instantaneously greatly exceed their mean 
value, resulting in spatially intermittent regions of high 
rates of turbulent dissipation within a turbulent flow field 
P, 0] ■ The statistical mean of the turbulent kinetic en- 
ergy dissipation rate equals the average flux of energy 
from the energy-containing large-scale eddies down to 
the smallest ones in case of statistically stationary tur- 
bulent fluid motion In the classical theory of tur- 
bulence, this small-scale end of the turbulent cascade 
is then prescribed by the Kolmogorov dissipation scale 
rjK which is derived as rjK = v'^^^ / {e)^^^ , where v is 
the kinematic viscosity of the fluid and (e) is the mean 
of the energy dissipation rate field which is given by 
e{x,t) = {i>/2){diUj + djUi)"^. But is rjK the smallest 
scale in a turbulent fiow? The Kolmogorov dissipation 
length contains the mean dissipation rate (e) and does 
not consider the strongly intermittent nature of e(a;, t). It 
seems therefore natural to capture for these fiuctuations 
in a dissipation scale definition and to refine the notion 
of one mean dissipation length to that of a whole contin- 
uum of local dissipation scales. The finest of those scales 
will then be associated with the steepest velocity gradi- 
ents, or alternatively, with the highest-amplitude shear 
and vorticity events in turbulence. This was done first 
within the multifractal formalism and demonstrated that 
local scales below rjK will exist 0, @, @i 01 ■ alterna- 
tive approach to such a continuum of dissipation scales 
was suggested by Yakhot lij . When connecting a lo- 
cal scale rj and the velocity increment across that scale, 
dr/U — \u(x + Tj) — u{x)\, the following relation 



(1) 



will be Rcri = rjSriu/v — 0(1). Such a Reynolds number 
is then also associated with fiuid motion at the crossover 
scales between the inertial cascade range and the viscous 
dissipation range. 

The probability density function (PDF) of the lo- 
cal dissipation scales, Q(?7), contains scales larger and 
smaller than the classical Kolmogorov dissipation length 
and can be calculated analytically on the basis of 
these ideas The analytic result was found re- 

cently to be in good agreement with a direct calcu- 
lation from high-resolution numerical simulation data 
of three-dimensional homogeneous isotropic box turbu- 
lence [13, 
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was obtained from the Navier-Stokes equations of fluid 
motion by a so-called point-splitting technique. Put into 
another perspective, a local Reynolds number obeying 
this scale and velocity as their characterisitic amplitudes 



The numerical simulations in [10] were 
conducted for relatively low Taylor-microscale Reynolds 
numbers (14 < Re\ < 151, where Rex = u'^X/v, 
A = y^2{15i/u'^) / (e) and the standard deviation of the 
turbulent velocity component is u'^) lacking even traces 
of inertial range. Still the scaling of moments of ve- 
locity gradients with Reynolds number agreed with the 
predictions which are obtained when using the scaling 
exponents of inertial range structure functions for high- 
Reynolds number flows. This somewhat unexpected re- 
sult suggested that the existence of a well-developed iner- 
tial range might not be essential for achieving an asymp- 
totic regime of the small-scale gradient statistics. 

In this Letter, we want to advance these ideas into two 
directions. We report the first experimental confirma- 
tion of the local dissipation scale concept in a laboratory 
fiow. Furthermore, we demonstrate the robustness of the 
concept beyond the idealized situation of homogeneous 
isotropic box turbulence. The flow at hand is a turbulent 
pipe flow in which a mean shear and thus an inhomoge- 
neous (radial) direction are present. Interestingly, we can 
also report that the monitored dissipation scale distribu- 
tions at the nearly isotropic pipe centerline and in the 
strongly anisotropic logarithmic layer are almost identi- 
cal for the range of Reynolds numbers that is spanned 
by the experiment. Additionally, they are in very good 
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agreement with those of the numerical simulations 11 1 
and the analytical result Our experiments span an 
Reynolds number range of 76 < Re\ < 214 . 

Description of the pipe flow experiment. The ex- 
periments were conducted using the Princeton Univer- 



TABLE I: Experiment conditions for measurements per- 
formed at y/R = 1. Ux is the time averaged streamwise 
velocity. 

Ren Rcl Re\ I7a;(m/s) it^m/s) I/(m) r7x(mm) (e)(m^/s^) 



sity/ONR Superpipe [I4I, which consists of a closed 


24000 


230 


76 


3.3 


0.12 


0.029 


0.47 


0.07 


return-pressure vessel containing a long test pipe down- 


28000 


290 


87 


3.9 


0.15 


0.029 


0.40 


0.13 


stream of flow-conditioning and heat-exchanging sec- 


35000 


340 


92 


4.9 


0.17 


0.029 


0.35 


0.21 


tions. The current set of measurements were performed 


44000 


440 


106 


6.0 


0.21 


0.030 


0.31 


0.39 


using air at room temperature and atmospheric pressure. 


52000 


560 


116 


7.2 


0.25 


0.033 


0.27 


0.62 


The test section used during these experiments consisted 


60000 


700 


124 


8.3 


0.29 


0.036 


0.25 


0.90 


of commercial steel pipe [13| with an average inner ra- 


70000 


780 


135 


9.5 


0.33 


0.036 


0.22 


1.31 



dius oi R = 64.9mm and an overall length of 200-D. Ex- 
periments were performed over an interval of pipe flow 
Reynolds numbers (based on pipe diameter D and area 
averag ed velocity Ub) of 2.4 x lO"* < i?eD = UbD/v < 
7.0 X 10'*. As detailed in [13,113], over this Reynolds 
number range, the flow at the test section of the pipe 
was fully developed, smooth walled turbulent pipe flow. 

Velocity measurements of streamwise velocity, Ux, 
were performed with a single sensor hot-wire probe 2.5/im 
in diameter and having a sensing length of I = 0.5mm. 
The wire was operated at an overheat ratio of 1.8. Ana- 
log output from the anemometer was low-pass filtered 
at half the sample frequency using an external filter be- 
fore being digitized using a 16-bit simultaneous sample 
and hold A/D board. Samphng rate was altered between 
lOOkHz and 200kHz in order to maximize the temporal 
resolution while not exceeding the measured frequency 
response of the probe. Data was sampled for 30 minutes 
continuously. During this period, the flow temperature 
was found to remain relatively constant, changing by less 
than rC. 

Hot-wire probe calibrations were performed at the pipe 
center line using a Pitot-probe/wall pressure tap combina- 
tion. Calibrations were performed before and after each 
30 minute measurement. Data sets where voltage drift 
was observed were discarded. 

Measurements were performed at two positions, the 
first located at the pipe centerline, and the second within 
the logarithmic layer at a distance from the wall of 
y = O.li?. The flow conditions for each measurement 
are provided in Tables [J and |TT] for the y/R = 1 and 
y/R = 0.1 measurements respectively. 

Confidence limits of the measurements. When resolv- 
ing the small scale structures of turbulence using hot-wire 
anemometry, the limitation imposed by the finite sensor 
length of the probe is an important consideration. In- 
sufficient spatial resolution will filter the energy at scales 
smaller than the probe. For the current measurements, 
the probe sensor length was of the order of the Kol- 
mogorov scale with a maximum of I « 477/f . Although 
it has been observed that the energy spectrum can be 
well resolved for turbulence scales up to 0.05Z [3], the 
impact of spatial filtering on the instantaneous velocity 
is not well understood, thus full confidence can only be 



TABLE II: Experiment conditions for measurements per- 
formed at y/R — 0.1. 

Ren Reh Rex Ux{m./s) it^m/s) -L(m) J7K(mm) (e)(m^/s^) 



35000 


1300 


155 


3.4 


0.41 


0.047 


0.19 


2.4 


52000 


1900 


185 


5.1 


0.59 


0.050 


0.15 


6.8 


70000 


2700 


214 


6.7 


0.76 


0.055 


0.13 


14 



given for measurements of 77 for scales larger than I. 

The scaled one-dimensional energy spectra, 4>(kx), are 
plotted in Fig. [T] for the centerline, y/R — 1.0, and in 
the logarithmic layer, y/R — 0.1. Streamwise wavenum- 
ber, kx, was estimated from frequency, /, using Taylor's 
frozen flow hypothesis through 2Tif /Ux- As can be ob- 
served in Fig. [1] the signal to noise ratio fell below one 
(as indicated by the inflection point in the measured spec- 
tra) at wavenumbers larger than the Kolmogorov range 
for all data sets. Note that the appearance of this noise is 
consistent with the so called noise typically observed 
in constant temperature anemometry. 

An isotropic estimate of the mean dissipation rate can 



TABLE III: Minimum values of r) which can be considered 
fully free of effects of 77; sensor spatial resolution and 77/ signal 
noise (defined using the wavenumber where the signal to noise 
ratio exceeded 10) . Scales are also shown normalized by rjo 



LRel'\ 


Rex 


Til (mm) 


Vi/Vo 


r?/(mm) 


Vf/Vo 


76 


0.5 


1.0 


1.4 


2.9 


87 


0.5 


1.2 


1.4 


3.5 


92 


0.5 


1.4 


1.5 


4.1 
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0.5 
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1.3 
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116 


0.5 


1.7 
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0.5 
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0.5 
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5.8 


155 


0.5 


2.2 
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4.5 
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0.5 


2.8 


0.7 


3.8 


214 


0.5 


3.2 


0.5 


3.0 
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FIG. 1: One-dimensional turbulence energy spectra measured 
for all cases. The legend indicates whether the data are ob- 
tained at the centerline y/R = 1.0 or in the logarithmic layer 
y/R^ 0.1. 
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FIG. 2: Measured dissipation spectrum for all cases, 
are the same as in Fig. [T] 
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be found through integration of the onc-dimcnsional dis- 
sipation spectrum 



(e) = 15i/ 



kl(j){kx) dkx 



(2) 



The measured dissipation spectra, shown in Fig. [2l indi- 
cate that the range of wavenumbers over which dissipa- 
tion occurred were well resolved during these measure- 
ments. Thus, the measured PDFs are expected to be 
largely insensitive to the influence of spatial filtering and 
instrumentation noise. The insensitivity of the results to 
changing Reynolds number and rjK in Fig. [3] confirm this 
expectation. Full confidence, however, can only be given 
to results for values of 77 larger than the scales at which 
these effects are expected to occur. These minimum con- 
fidence limits are quantified in Table UTIl 



FIG. 3: Measured PDFs of the local dissipation scales for 
all Reynolds numbers. Results from the logarithmic layer 
y/R ~ 0.1 are shifted vertically by a decade for clarity. 



Distribution of local dissipation scales. The probability 
density function (PDF) of 77, Q{i]) can then be calculated 
by a saddle-point approximation ^ within the range < 
1] < L as 



Qiv) = 



dx 



7777(6 log (L^))^/^ J-c 

(log (^V2xReL {v/Lf 



X exp 



Ah\og{L/vi) 



(3) 



where a = 0.383, h — 0.0166 11], L is the integral 
length scale of the turbulence and Re^ = {\ux{x -I- L) 
)^l'^Llv as used in 



{x)Y 



Uj. The PDF was calcu- 



lated from each velocity time series using the follow- 
ing procedure. First, the velocity difference Au(At) = 
|J7a;(i-l- Ai) — C/a;(t)| was calculated at each discrete time, 
t. Next, the value Au(Ai)C/a;Ai/z/ was calculated and 
the instances where the result was between 0.9 and 2 
were counted as occurrences of dissipation at a scale 
77 = Ate/ 2:. This process was repeated for all Ai up to 
^tU X = Finally Q{r\) was found from the count of 
occurrences by normalizing such that / Q(r\)dr\ — 1. 

The results are shown in Fig. [31 normalized by the 
value 770 



LRe 



-0.72 



ll| . which is very close to r\K. 
At each measurement location there was excellent col- 
lapse of Qi'q/'qo) for all measured Reynolds numbers. 
This Reynolds number independence indicates that the 
high-wavenumber end of the universal equilibrium range 
could establish universality at lower Reynolds numbers 
than previously expected. This appears to be true even 
though the scale separation between the energy contain- 
ing eddies and dissipation scales is still relatively small 
and before even the establishment of an inertial subrange 
(see Fig. [J). 
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FIG. 4: Comparison between Qirj/rjo) measured at pipe cen- 
terline, y/R = 1.0, measured at y/R = 0.1, numerical results 
of [ll|| and theoretical distribution of [^. 



Direct comparison is made in Fig. |4] between Q^rj/ijo) 
measured in the nearly isotropic turbulence at the pipe 
centerline and the anisotropic log-layer. Here, one can 
see nearly identical distributions at all but the smallest 
values of 77/7/0 where a slightly increased probability was 
observed at y/R = 0.1. Such a good agreement indicates 
that the structure of turbulence at the largest scales has 
little effect on the organization of the smallest scales. 

Comparison is also made in Fig. |4]between the present 
experimental results, the numerical results of [Tlj and the 
theoretical distribution of [1] . There is a very good agree- 
ment between the numerical and experimental results (at 
the center line), also in the tails of the PDF at large 77. 
The deviations from the theoretical distribution on both 
sides can be due to the limits of a saddle-point approxi- 
mation for the evaluation of the Mellin transform in 
The agreement between the numerical and experimental 
results supports the concept of the universality of the 
distribution of the dissipation scales. 

As mentioned above. Fig. [4] displays a further prop- 
erty. The left tail of the PDF Qirj/rjo) becomes slightly 
fatter for the logarithmic layer in comparison to that of 
the center line. Since we found that the data for dif- 
ferent Reynolds numbers collapse well (see Fig. [3]), this 
property could be interpreted as a statistical fingerprint 
of bursting structures. So-called packets of hairpins have 
been observed recently in the logarithmic layer of other 
wall flows [l3|. They were found to increase the level 
of small-scale intermittency which would be in accor- 
dance with the slightly smaller measured local dissipa- 
tion scales. However, a clear disentanglement requires 
two things, further data records at higher Re combined 
with experiments in different flows. 

Conclusions. The probability density function of the 
local dissipation scales was measured in a turbulent pipe 



flow over a range of Reynolds numbers at the pipe cen- 
terline and in the logarithmic layer. For the first time, 
we measured this probability density function experi- 
mentally in a turbulent shear flow. Our results indicate 
that the distribution is basically independent of both the 
Reynolds number and the degree of anisotropy of the 
large scales of turbulence. We find a very good agreement 
with the numerical simulations and the theory. This ro- 
bustness of the results with respect to flow-type (or de- 
gree of anisotropy) suggests a universality of the smallest- 
scale fluctuations in turbulence. One could thus indeed 
conclude that the turbulent dynamics at the finest scales 
is already in an asymptotic state although the Reynolds 
numbers remain moderate and an inertial cascade range 
is absent or very small. It is well-known that to observe 
inertial range statistics in wall-bounded flows requires 
large Reynolds numbers [l^ . This fact adds further value 
to the present study which establishes a connection be- 
tween shear flow and isotropic turbulence at their small 
scales. 
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